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Abstract. We report the discovery of the gamma-ray burst 
GRB 000131 and its optical afterglow. The optical identifica- 
tion was made with the VLT 84 hours after the burst follow- 
ing a BATSE detection and an Inter Planetary Network local- 
ization. GRB 000131 was a bright, long-duration GRB, with 
an apparent precursor signal 62 s prior to trigger. The after- 
glow was detected in ESO VLT, NTT, and DK1.54m follow- 
up observations. Broad-band and spectroscopic observations of 
the spectral energy distribution reveals a sharp break at optical 
wavelengths which is interpreted as a Lya absorption edge at 
6700 A. This places GRB 00013 1 at a redshift of 4.500±0.015. 
The inferred isotropic energy release in gamma rays alone was 
~ 10 54 erg (depending on the assumed cosmology). The rapid 
power-law decay of the afterglow (index a=2.25, similar to 
bursts with a prior break in the lightcurve), however, indicates 
collimated outflow, which relaxes the energy requirements by 
a factor of < 200. The afterglow of GRB 000131 is the first to 
be identified with an 8-m class telescope. 

Key words: cosmology: observations - gamma rays: bursts - 
ISM: dust, extinction 



1. Introduction 
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Since the first detections of soft X-ray (Costa et al., 1997), op- 
tical ( van Paradijs et al., 1997] ), and radio (Frail and Kulkarni, 

1997) afterglows to gamma-ray bursts (GRBs), 19 optical af- 
terglows have been detected. 

For the current sample of 12 GRBs with well-determined 
redshifts, the median redshift is z me( jian =1-1 with a very large 
root-mean-square variation of ~0.8. While a photometric red- 
shift of w 5 has been proposed for GRB 980329 (Fruchter, 
1999), the highest spectroscopic redshift determination for a 
GRB so far is z = 3.418 for GRB 971214 (Kulkarni et al., 

1998) . This redshift was determined from the spectrum of a 
faint Lyman Break galaxy that was found at the position of the 
optical afterglow (see also Odewahn et al., 1998). The high- 
est redshift determined directly from absorption lines in the 
spectrum of an optical afterglow is z = 2.0404 ± 0.0008 for 
GRB 000301C dSmette et al., 20001; llensen et al., 2000b. Sev- 



eral authors have speculated about the possibility of detecting 
GRBs at even higher redshifts than achieved so far. Wijers et 
al. ( |1998h and Lamb & Reichart ( |2000| ) estimate that GRBs 
should be detectable at very high redshifts (z > 5) based on the 
assumption that GRBs are associated with star formation and 
by using models for the cosmic star formation rate as a function 



of redshift. Blain & Natarajan ( 2000 ) propose to use the oppo 



site strategy, namely to use the observed redshift distribution of 
GRBs to measure the (uncertain) cosmic star formation rate as 
a function of redshift. 
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Table 1. Spectral data from BATSE 



Peak Flux: 7.89 ± 0.08 phcirT^ s _1 (46-313 keV; 1.024 s) 

Fluence: (3.51 ± 0.08) Xl0~ 5 erg cm" 2 (26-1800 keV; 137 s) 

Duration: T90 = 96.3 s (error ~ 1 s) 
T50 = 30.7 s (error ~ 1 s) 

Time-averaged E pcak = 163 ± 13 keV 

spectral Low-energy index — —1.2 ±0.1 

parameters: High-energy index — —2.4 ±0.1 



Unfortunately, afterglows have been detected in only < 40 
percent of all well-localized (SAX, RXTE, IPN) GRBs (e.g., 
Fynbo et al. 2000). It is not clear why most GRBs with no ap- 
parent optical afterglow ('dark bursts') are so abundant. Possi- 
ble explanations are rapid decay (3root et al., 1998), reddening 
in the immediate environment and host galaxy (eg., Jensen et 
al. (2000)), intrinsic optical faintness ( Taylor et al., 2000), very 
high redshifts ( fr^amb & Reichart, 2000 ) or a combination of 
these effects. To shed light on dark bursts, one obvious strategy 
is to set deeper optical limits by using large telescopes and to 
determine the decay slopes and colours (see Fynbo et al. 2000). 

In this Letter we present the first attempt to identify the af- 
terglow of a GRB with an 8-m class telescope. The attempt was 
successful in that the observations led to the discovery of the 
afterglow and the subsequent determination of the decay slope, 
spectral energy distribution and record high redshift. In the fol- 
lowing we document the gamma-ray, optical, and infrared ob- 
servations, and discuss the properties of the gamma-ray burst 
and afterglow. 
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Fig. 1. Gamma-ray light curve of GRB 000131 recorded by the 
BATSE Large Area Detectors (LADs). b shows continuously 
recorded data with 1024 ms time resolution, a shows an en- 
largement of the most intense interval using triggered data with 
64 ms time resolution. In both cases, data from the two most 
brightly illuminated LADs have been summed. Energy chan- 
nels for the different data types and detectors have been chosen 
to approximately match (within a few percent) the energy range 
60-320 keV. 



2. Detection of the gamma-ray burst 

GRB 000131 was observed by Ulysses, KonusAVind, NEAR, 
and CGRO-BATSE on 2000 Jan. 31.624 UT. After 56 hours, it 
was localized via Interplanetary Network (IPN) timing to two 
alternate 55 sq. arcmin error boxes (Hurley et al., 2000a), one 
of which was ruled out by the independent BATSE location 
(Kippen, 2000). From the Ulysses, KonusAVind, and NEAR 
data its 25-100 keV fluence was ~ 10~ 5 erg/cm 2 . 

BATSE detected the event (Trigger #7975) in a partial data 
gap, so the standard catalog data products such as flux, fluence, 
and duration are not available. However, by analyzing other 
data types from BATSE the relevant parameters could be esti- 
mated. The results are given in Table [j]. The peak flux is in the 
top 5%, and the fluence in the top 7% of all BATSE GRBs. The 
spectrum is well fit by a standard GRB model (Band et al., 
1993 ) and shows a typical spectral evolution (P reece et al., 
2000). 

The BATSE light curve is shown in Fig. [j]. The event con- 
sists of several spikes, most of which show a clear asymmetry, 
the leading edge being the steepest. A weak pulse lasting « 7 s 
is observed 62 s prior to the burst trigger. Its arrival direction 
is consistent (within a 12 degree error) with that of the main 
pulse. This spike is therefore most likely part of the burst. In 
this case the overall duration is about 170 s, i.e. significantly 
longer than measured by T90. The presence of this 'precursor' 



is reminiscent of that displayed by GRB 991216 (Kippen et al., 
1999; hurley, 2000a| ). However, in the case of GRB 000131, no 
long-lasting tail (7-ray afterglow) is observed. 



3. Localization and identification of the afterglow 

The late localization of the GRB complicated the detection of 
the optical afterglow, expected to be significantly fainter than 
the limit of existing sky surveys. Furthermore, the (presumed 
power-law) decay of the afterglow would be so slow that obser- 
vations separated by several days would have to be compared 
before any probable candidate source could be established. 

Images of the Inter Planetary Network (IPN) error box were 
obtained with the FORS 1 instrument on Antu (ESO VLT UT1), 
starting 84 h after the burst. Four 120 s exposures in the B,V 
and R bands of each of two FORS 1 fields, covering the error 
box (see Fig. ^), were acquired under good seeing conditions 
(see Table [| for a log of the observations). On the same night 
we acquired I-band exposures with a total integration time of 
3600 s, using DFOSC on the Danish 1.54-m telescope on La 
Silla. From this first set of exposures no candidate optical tran- 
sient could be identified. A subsequent set of images with the 
same exposure times was acquired 135 h after the burst, under 
less favorable seeing conditions. Comparing these two epochs, 
one source located at R.A. = 6 h 13 m 31!08, Dec. = -51°56'4T/7 
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Table 2. Journal of our observations of GRB 000131. 




GRB 000131 
(06;1.3:31.8,- 



-51:56:4-1 



Fig. 2. The rhombus shaped IPN error box with the two square 
FORS1 fields superimposed. The location of the OT (not 
present in this Digital Sky Survey image) is marked. 

(J2000) was found to have declined by about 1 . 1 mag in the 
R-band. In the first epoch images this source was also de- 
tected in the V-band, but not in the B-band. Its proximity (1 .'1) 
to the center of the error box, in combination with a fading 
of about half a magnitude per day, which is typical for GRB 
afterglows after some days, made this obje ct the likely after- 
glow of GRB 000131 ( [Pedersen et al„ 2000| ). The afterglow of 
GRB 000131 was t he second to be identified based solely on 
an IPN localization (|Hurley, 2000^). 



4. Follow-up observations and data analysis 

Additional images of the field were obtained on 2000 Feb. 8.1 
UT. The source was detected in a deep FORS 1 R-band expo- 
sure and was found to have faded by 1.9±0.2 mag relative to 
the first epoch confirming its transient nature. As the source, 
based on the V, R and I photometry, appeared to be very red, 
we also acquired near-infrared (JHKs) images with the ESO 
New Technology Telescope (NTT). A final set of R and I band 
images were acquired with FORS1 on 2000 Mar. 5.0 UT. A 
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Table 3. Magnitudes of internal reference stars. 



Star 


B 


V 


R 


I 


J 


H 


K 


A 


18.97 


18.49 


18.13 


17.77 


17.38 


17.01 


16.88 


B 


21.60 


20.07 


18.86 


17.48 


16.12 


15.53 


15.25 


C 


20.79 


19.90 


19.32 


18.86 


18.19 


17.59 


17.53 


D 


21.15 


19.68 


18.65 


17.67 


16.54 


15.92 


15.73 



log of our photometric observations and the standard calibrated 
magnitudes of the afterglow is given in Table ||. All image re- 
ductions were performed in the IRAF[] environment. 



4.1. Image analysis and photometry 

The afterglow was detected at low signal-to-noise in most of 
the images, which made it essential to use PSF photometry to 
derive the magnitudes. We used DAOPHOT-II ( fctetson, 1987b 
within IRAF to derive PSF magnitudes of all point sources in 
the images. In general, when deriving PSF magnitudes, the cen- 
ter and amplitude (magnitude) of the PSF is fitted. For objects 
close to the detection limit, fitting the centroid introduces a bias 
in the magnitude of up to one magnitude, because of a tendency 
to fit the PSF to a nearby noise spike. If instead DAOPHOT is 
used on objects with accurately known pixel coordinates, the 
magnitude can be derived without re-centering the object dur- 
ing the PSF fit. In addition to removing the magnitude bias, de- 
riving the PSF magnitude without re-centering also has the ad- 
vantage that magnitudes down to the 2-er level can be obtained. 
We could therefore use the accurate position of the afterglow, as 
determined from the first epoch images, to derive reliable mag- 
nitudes from the observations at later epochs. Limiting magni- 
tudes were obtaind by adding a set of artificial point sources 
to the images at locations of apparently blank sky, setting the 
limiting magnitude at the level where DAOPHOT ALLSTAR 
could recover 95% of the point sources. 



IRAF is the Image Analysis and Reduction Facility made avail- 
able to the astronomical community by the National Optical Astron- 
omy Observatories, which are operated by AURA, Inc., under contract 
with the U.S. National Science Foundation. 
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Fig. 3. R-band image of the GRB 000131 afterglow from 
Feb. 4. The afterglow is marked with two bars and the four 
comparison stars are labeled according to Table |[ North is up 
and East is to the left. 



4.2. Optical images 

The optical observations were all acquired on nights with pho- 
tometric sky quality. The photometry was calibrated to standard 
star fields (Landolt, 1992) through CCD aperture photometry 
relative to the comparison stars marked in Fig. || and listed in 
Table [3} The errors in the magnitudes of the reference stars are 
dominated by the determination of the photometric zero points. 
Extinction correction was applied using observatory standard 
extinction values. Because the standard star exposures were ac- 
quired at approximately the same airmass as the object expo- 
sures, corrections were in all cases smaller than the individual 
photometric error of the afterglow magnitude. From the errors 
in the color transformation, it is estimated that the photometric 
zero point errors are smaller than 0.02 mag. In the R and I im- 
ages from Mar. 5, there is no detection of a source within 2" of 
the afterglow location. This implies an upper limit to the host 
galaxy magnitude of R = 25.7 and I = 24.8. 



The frames were reduced by first subtracting a mean sky, 
obtained from frames acquired just before and after the source 
frame. Before using the frames for sky subtraction, stars were 
eliminated with an automatic star finder coupled with a back- 
ground interpolation algorithm. Then, a differential dome flat- 
field correction was applied, and the frames were registered to 
fractional pixels and combined. We calibrated the photometry 



with standard stars taken from Persson et al. (1998), acquired 
directly before and after the source observations. The standards 
were placed in five different positions on the detector, and re- 
duced in the same way as the source frames. Standard star and 
source photometry was corrected for atmospheric extinction 
using the me an ESQ La Silla extinction coefficients given in 
Engels et al. ( 1981 ). Formal photometric accuracy, as judged 
from the standard deviation of the standard star observations, 
is about 0.03 mag. 

At the time of observations, about 8 days after the IPN de- 
tection, the afterglow was very faint and thus only marginally 
detected in H (3 a) and K (2.5 a). In J there is apparently a 
source, but it is displaced about 1 pixel from the location of the 
afterglow. A PSF magnitude of this source could only be de- 
rived if re-centering during the PSF fit was allowed. To main- 
tain consistency in our magnitudes, we have therefore chosen 
to disregard this apparent source and only give the proper limit- 
ing magnitude. The H and K magnitudes are both very close to 
the detection limit. It was verified that magnitudes at this level 
can be reliably derived by performing photometry of a large 
number of artificially added stars of the same magnitude. 

4.4. Spectroscopy 

Starting 2000 Feb. 8.09 UT, we obtained spectroscopic obser- 
vations, using FORS 1 configured with a 300 lines/mm grism, 
a cut-on order sorter filter and a 0V7 slit. The wavelength range 
covered was 3800 A - 8200 A, with a resolution of 9 A and 
with order overlap from about 7600 A. The total integration 
time was 10800 s, divided into 6 exposures of 1800 s. During 
the spectroscopic observations, the seeing ranged from 0'.'55 to 
079. The combined spectrum was flux calibrated relative to the 
flux standard LTT3218. The continuum is not well detected at 
visible wavelengths but is weakly detected in the red part of 
the spectrum, with a signal-to-noise around 1 per spectral bin 
of 2.5A in regions that are not dominated by sky lines. The 
spectral slope is consistent with the broad-band photometry. 
The spectrum does not show signs of statistically significant 
emission lines, as could be expected from a host galaxy. 



4.3. Near-infrared images 



5. Discussion 



Near-infrared (IR) J, H, and Ks-band images were obtained at 
La Silla with the NTT and SOFI on 2000 Feb. 8.1 UT. SOFI 
is equipped with a Hawaii 1024x 1024 pixel HgCdTe detector, 
and we used a plate scale of 0'.'29 which gives a field-of-view 
of roughly 4.'9x4f9. Each image comprises several tens of el- 
ementary co-added frames acquired by randomly dithering the 
telescope several arcsec once a minute. 



The temporal and spectral behavior of optical gamma-ray burst 
afterglows are in general observed to follow a power law, such 



that F u (t,v) oc r 



(a, j3 > 0), where t is the time 



elapsed since the GRB event. The fireball model provides the 
theoretical framework for this behavior and also gives the re- 
lation between the t emporal slope a and the spectr al slope (3 
in various regimes ( Sari et al., 1998 ; Piran, 1999 ; Meszaros 
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Table 4. GRB afterglows with rapidly decaying light curves. 



et al., 1998). If the afterglow is collimated (a jet), the temporal 
slope is modified and may change during the transition from 
the observer being inside the jet cone to bein g outside the je t 
cone, caused by the slowing-down of the jet ( Rhoads, 1999 ). 
This results in a break i n the light curve, su ch as those ob- 
served in GRB 980519 ([launsen et al., 2000|), GRB 990123, 



GRB 


a 


P 


tbrcak 


R(t=4 d) 




980326 


2.10 ±0.13 


0.66 ±0.70 




26.4 


20 - 


980519 


2.22 ±0.04 


0.80 ±0.08 


0.55 ±0.17 


24.7 




990510 


2.18 ±0.05 


0.61 ±0.12 


0.80 ±0.35 


22.1 




991208 


2.20 ±0.20 


0.75 ±0.03 




20.2 


22 - 


000131 


2.25 ±0.19 






23.6 


a; 


00030 1C 


2.29 ±1.00 


0.55 ±0.04 


4.39 ±1.52 


20.6 


24 - 



GRB 990510 (Stanek et al., 1999 



GRB 990705 flMasetti et al., 2000a|), GRB 991216 (Halpern 
et al., 2000), and GRB 0001 3 1C Jensen et al., 2006) Masetti 



Harrison et al., 1999), 



et al., 2000b). In the following, we will base our discussion on 
the assumption that the afterglow can be described within the 
framework of the fireball model. 
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Fig. 4. R-band light curves of GRB afterglows with steep late- 
time decay. The solid line is a power-law fit to the GRB 
000131 data (filled circles). The non-detection on March 5 is 
marked as an upper limit at log £=1.52. Other afterglows are 
(from left to right) GRB 980326, GRB 980519, GRB 990510, 
GRB 00030 1C and GRB 991208 (which is shifted 0.4 in 
log(days)). Pre-break data points are omitted. 



5.1. The light curve of the afterglow 

The three epochs of R-band photometry can be fitted by a 
power-law decline with a = 2.25 ± 0.19 and a \ 2 = 0.06, as 
shown in Fig. |j. A power-law decline with this slope is typical 
for post-break evolution of GRB afterglow emission. A value 
of 2.0 < a < 2.5 has been observed in several other after- 



glows, namely GRB 980326 ( proot et al., 1998|), GRB 9 80519 
dJaunsen et al., 2000|), GRB 9905 10 ( |Stanek et al., 1999| ; Harri- 
son et al., 1999; |Holland, 2000|) GRB 991208 (Hurley et al., 
1999; |Castro-Tirado et al., 2000| ), and GRB 000301C (e.g., 
Jensen et al. 2000). The light curve of GRB 000131 is plot- 
ted together with the light curves of these afterglows in Fig. [5J. 
The main characteristics of these systems are summarized in 
Table Q In all cases the favoured interpretation of such a steep 
slope is post-break decay of a collimated outflow. The light- 
curve breaks predicted in this scenario have indeed been ob- 
served (see Holland et al. 2000) and are believed to have oc- 
curred prior to the first observations in the other systems with 
a > 2. The simplest interpretation of the GRB 000131 after- 
glow light curve is therefore that of a collimated outflow seen 
after the jet has slowed down. 



5.2. Spectral energy distribution 

The optical afterglow was detected in VRIHK, but at differ- 
ent epochs. By assuming achromatic evolution of the afterglow 
(which is consistent with other afterglows), following the a = 
2.25 power-law decline, the spectral energy distribution can be 
calculated for a given epoch. In Fig. |5] we plot the derived spec- 
tral energy distribution at 3.50 days after the burst trigger, i.e., 
at the time of the first VRJ observations. Also given are the up- 
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Fig. 5. The spectral energy distribution, as derived from broad 
band photometry. The errors of the H and K fluxes include the 
formal error from the extrapolation of the light curve back to 
£=3.5 days. A fit by a power law spectrum with Lyman forest 
absorption and SMC reddening is shown as a dashed line. This 
yields Ay = 0.18, when an intrinsic spectral slope, = 0.70, 
and a redshift of 4.5 is assumed. The solid line shows the cor- 
responding spectrum with its Lyman absorption edges. 



per limits derived from our B and J observations. The spectral 
energy distribution has been corrected for Galactic reddening, 
using E (B _ V )=0.056 ( |Schlegel et al., 1998| ). 



For steep-decay afterglows, like those shown in Fig. the 
large values of a and the relatively small values of (3 favour a 
scenario involving a sideways expanding jet. In this case spec- 
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tral slopes between (a — 1) /2 and a/2 are expected in the fire- 
ball model, depending on the value of the cooling frequency. 
With a = 2.25, this implies 0.63 < (3 < 1.13, with a prefer- 
ence for the low value (cf. Table Q). 

The spectral energy distribution shown in Fig. || does not 
resemble a power-law with a single index (3. This is not an arti- 
fact of extrapolating the H and K band data points from Feb. 8 
to Feb. 4, as a power-law fit to the V, R, and I data points only 
results in a x 2 = 7.1 and (3 = 4.90, which is much steeper than 
the value of about 1 typically observed in GRB afterglows (see 
Table Q). Moreover, for the H and K fluxes to be in accordance 
with this fit, an unphysical value of a of about 7.0 in the near- 
IR is required. Thus, a power-law spectral energy distribution 
is ruled out. 

We have explored whether the strongly curved shape of the 
spectral energy distribution can be explained by reddening in 
the host galaxy (as in eg. GRB 000301C, Jensen et al. 2000, 
or in GRB 971214, Dal Fiume et al. 2000). We find that no 
physically plausible reddening laws can transform a power-law 
spectrum into the observed shape. 

The most likely interpretation of the spectral break is there- 
fore the onset of Lyman forest blanketing, hence implying that 
GRB 000131 was at a high redshift (z>4). To examine this in- 
terpretation further we have fitted the spectral energy distribu- 
tion with power-laws modified by the effects of Lyman forest 
blanketing and internal reddening in the host galaxy. This tech- 
nique was first used by Fruchter (|1999| ) for GRB 980329 and 
discussed in detail by Reichart ( 2000| ). For the reddening we 
use the SMC extinction law assumed to best represent a chem- 
ically less-evolved environment at redshifts of 4 to 5. To model 
the effects of Lyman forest blanketing we follow M0ller and 
Jakobsen (1990): for each pair of values of (3 and z, the visual 
rest-frame absorption Ay and the \ 2 of the fit is obtained. We 
find that all values of (3 allowed by the fireball model are possi- 
ble, with a preference for low values. The corresponding range 
of internal reddening is 0.11 < Ay < 0.20, with a preference 
for high values. For values of around 0.65, the possible range 
(2 a) of redshifts is 4.38 < z < 4.70, while the redshift is 
4.55 < z < 4.69 for (3= 1.0. This result is not affected by the 
choice of reddening law. 

5.3. The spectroscopic redshift of GRB 000131 

Our analysis of the photometric observations implies the pres- 
ence of a Lya absorption edge in the range 6500 A to 6900 A. 
This implication is confirmed by the spectroscopic observa- 
tions. After smoothing the FORS1 spectrum to a resolution of 
30 A a clear indication of such an absorption edge at 6700 A 
is revealed. The recorded spectrum was very faint, at the level 
of less than 2% of the continuum of the night sky, which is 
dominated by emission lines through the red part of the spec- 
trum. The noise in the spectrum is therefore dominated by un- 
certainty in the sky subtraction. An improved representation 
of the continuum was obtained by rebinning the spectrum to 
a resolution of 200 A, omitting spectral bins coincident with 
sky lines. The smoothed and the rebinned spectra are shown in 
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Fig. 6. The spectrum of the GRB 000131 afterglow: a 
smoothed to a resolution of 30 A, b rebinned to a resolution 
of 200 A. The dotted curve shows the noise per bin, while the 
dot-dashed line is the model spectrum shown in Fig. ||J 



Fig. |6] together with a model of the spectrum, assuming a red- 
shift of 4.5, (3 = 0.70 and a reddening, Ay = 0.18, as derived 
from a fit to the photometry (see Fig. ^|). 

The model spectrum is normalized to the R-magnitude of 
the afterglow, obtained immediately before the spectroscopic 
observations began, and corrected for slit losses of 35%, as 
would result from using a ff.'l slit in a seeing of about 0'.'7. 
The model spectrum is seen to be in very good agreement with 
the binned spectrum for wavelengths below 7200 A. The av- 
erage observed flux in the Lyman forest region (between Lya 
and Ly/3) is 55 ± 15 nJy, which is in reasonable agreement with 
the level of 79 nJy predicted by the model. Beyond 7200 A the 
spectrum is dominated by an atmospheric band and strong sky 
lines, which are not well resolved at a spectral resolution of 
9 A. This effectively renders the spectrum useless in this spec- 
tral region. 

As the spectral region from 6590 A to 6810 A is essen- 
tially free from sky lines, the location of the absorption edge 
is well defined (see Fig. ^). We measure its location to be at 
6701 ± 2.5 A from which a redshift of 4.511 ± 0.002 is in- 
ferred. The reality of the absorption edge and the interpretation 
that it is due to Lya is most convincingly seen by integrating 
the flux of the spectrum and the model, as shown in Fig. [7| 
The integrated spectrum is very smooth and follows the model 
nicely on the red side of the edge, while there are significant 
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deviations on the blue side of the edge. Comparing 100 A in- 
tervals on the blue and red side of the absorption edge, the stan- 
dard deviation of the spectrum is found to be 65% larger on the 
blue side, consistent with the interpretation that the blue side is 
located in the Lyman forest. Hence, the VLT spectrum is in ex- 
cellent agreement with the photometric observations and pro- 
vides independent evidence that the redshift of GRB 000131 is 
z w 4.5. 



In the VLT spectrum of GRB 000301C, Jensen et al. (gOOfj) 
detected a strong damped Lya absorption line at the redshift of 
the GRB. If a similar damped Lya absorption line is present in 
the spectrum of GRB 000131 the redshift as determined from 
the spectral break will be slightly overestimated since the red 
wing of the damping profile will move the spectral break 10- 
20 A (depending on the H I column density) towards the red. 
Therefore, the inferred redshift from the spectral break depends 
on the assumed H I column density of the GRB self absorp- 
tion. With the signal-to-noise of the spectrum, it is not possi- 
ble to constrain the H I column density. An indirect indication 
that GRB 000131 had significant self absorption is obtained 
from the estimate of Ay - As suming an SMC extinction law and 
minimum dust-to-gas ratio ( Pei, 1992 ) and the maximum value 
of Ay, as derived from the fit of the spectral energy distribu- 
tion, an H I column density of up to 10 22 cm -2 is estimated. 
An estimate of the likely error of the redshift may therefore 
be obtained by fitting a H I line corresponding to this column 
density to the spectrum. We derive a redshift of 4.490 ±0.002 
for this column density, which allows us to conservatively con- 
clude that the redshift of GRB 000131 is z = 4.500 ± 0.015. 

6. Conclusions 



Based on the temporal and positional coincidence of the tran- 
sient object with GRB 000131, and the fact that it was lo- 
cated at a cosmological distance, we conclude that this tran- 
sient object was the optical afterglow of GRB 000131. Assum- 
ing a Hubble constant of 65 km s -1 Mpc -1 , fi = 0.3 and 
£!a = 0.7, a redshift of 4.50 implies a distance modulus of 
48.24. Hence, the total energy release assuming isotropic emis- 
sion is k 1.1 x 10 54 erg. This is only a factor of ~3 smaller than 
the isotropic equivalent energy for the most energetic event so 



far, GRB 990123 (Andersen et al., 1999; Kulkarni et al., 1999) 



However, the combination of the decay slope a and the con- 
straints on the spectral slope j3 indicate that the afterglow was 
due to a collimated jet. The upper limit on the break epoch, 
tbreak < 3.5 days, implies a lower limit on the jet opening an- 
gle of 8 > 7° n 1 / 8 > 7°, where n is the density of the ambi- 
ent medium (in units of cm -3 ), which is larger than 1 in star- 
forming regions. From the lower limit on the opening angle we 
infer a lower limit on the released energy of ~ 5 x 10 51 erg. 

This work shows that 8-m class telescopes may be used 
successfully to detect faint R > 23 optical afterglows even in 
fields with very bright stars. This is necessary in order to re- 
solve why > 60% of all attempts to detect optical after glows 
of GRBs are currently unsuccessful ( Fynbo et al., 200C| ). The 
apparent brightness of the GRB 000131 afterglow was similar 
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Fig. 7. a A region of the spectrum, centered on the Lya absorp- 
tion edge and smoothed to a resolution of 13 A. The dotted line 
gives the noise per bin, while the dashed curve is a model ab- 
sorption edge spectrum corresponding to a redshift of 4.511, 
and redshifts of 4.509 and 4.513 (dot-dashed curves), b The 
integrated flux of the full resolution spectrum and the model 
absorption edge spectra shown in a. The integrated flux is nor- 
malized to a common level at 6750 A . 



to that of GRB 990123. GRB 000131 does therefore provide 
the first observational evidence that it is possible to obtain high 
resolution optical spectra of GRBs at very high redshift, if the 
afterglow is identified at an early time. High redshift GRBs will 
therefore no doubt prove to be an extremely valuable tool for 
not only the understanding of the GRB environment but also 
for the study of the Lyman forest and of cosmology. 
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